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Abstract 



It is often useful in multivariate time series analysis to determine statistical causal 
relations between different time series. Granger causality is a fundamental measure 
for this purpose. Yet the traditional pairwise approach to Granger causality analysis 
may not clearly distinguish between direct causal influences from one time series to 
another and indirect ones acting through a third time series. In order to differentiate 
direct from indirect Granger causality, a conditional Granger causality measure in 
the frequency domain is derived based on a partition matrix technique. Simulations 
and an application to neural field potential time series are demonstrated to validate 
the method. 
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1 Introduction 



The concept of causality introduced by Wiener (Wiener, 1956) and formulated 
by Granger (Granger, 1969) has played a considerable role in investigating 
the relations among stationary time series. The original definition of Granger 
(1969), which is well named as Granger causality, refers to the improvement in 
predictability of a series that derives from incorporation of the past of a second 
series, above the predictability based solely on the past of the first series. This 
definition only involves the relation between two time series. As pointed out 
by Granger (Granger, 1969; 1980), if a third series is taken into account, a 
spurious or indirect causality due to the third series may be detected. Then he 
defined a prima facie cause (Granger, 1980): Y is said to be a prima facie cause 
of X if the observations of Y up to time t (Y(t) : r <t) help one predict X(t+ 
1) when the corresponding observations of X and Z are available (X(r), Z{r) : 
t < t). We refer to this idea as conditional Granger causality since it gives a 
measure of causality between two time series, X and Y, conditional on a third, 
Z. Evaluation of this conditional Granger causality in the time domain is fairly 
straightforward through comparison of two predictions of X(t + 1), one when 
(X(t), Z(t) : t < t) are given, the other when (X(r), Y(r), Z(r) : t < t) are 
given. However, evaluating causality by frequency decomposition may allow 
more meaningful interpretations in cases where oscillations are involved. 

After giving clear measurements of linear dependence and feedback between 
two blocks of time series (Geweke, 1982), Geweke also presented a measure of 
conditional linear dependence and feedback (Geweke, 1984). Both a time do- 
main measure, consistent with that of Granger, and its frequency decomposi- 
tion were given. Although Hosoya presented some improvements on Geweke's 
methods (both bivariate (Hosoya, 1991) and conditional versions (Hosoya, 
2001)), they have not been widely accepted because his time domain imple- 
mentation departs from Granger's original idea, and its physical interpretation 
is less clear. 

We point out that Geweke's use of the term "feedback" is equivalent to "causal- 
ity" in the present discussion. In applying Geweke's frequency-domain condi- 
tional Granger causality measure to neural field potential data, we have found 
that negative values, which have no meaning in terms of causality, may occur 
at some frequencies. This finding casts doubt on the applicability of Geweke's 
method for neural time series analysis. We believe that the negative values re- 
sult from the lack of identity of estimates of the same spectrum when different 
autoregressive (AR) models are used. This non-identity of different estimates 
of the same spectrum is a general practical problem in numerical analysis that 
causes errors in Geweke's implementation because it requires the estimates to 
be identical. In this paper, we employ a partition matrix method to overcome 
this problem. Comparison of the results from our procedure with Geweke's 
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original procedure, clearly shows the validity of the current procedure. In the 
following sections: we first provide an introduction to Granger causality; then 
present an overview of Geweke's procedure on conditional causality, pointing 
out the importance of obtaining a correct measure; and then derive our pro- 
cedure. Finally, results of simulations and application to neural field potential 
time series data are provided. 



2 Background 

Consider a multiple stationary time series of dimension n, W = {w t }. The 
series has the following moving average representation with use of the lag 
operator L: 



where E(e t ) = 0, var(e t ) = £ and A = I n , the nxn identity matrix. Assume 
there exists the autoregressive representation: 



where B = I„. 

Suppose that w t has been decomposed into two vectors x t and y t with k 
and / dimensions respectively: w t = (x£, y' t )', where the prime denotes matrix 
transposition. Denote W t _i as the subspace generated by {w s ; s < t — 1}. 
Define S x = var(x t |X t _ 1 ), E 2 = var(x t |X t _ 1 , Y t _0, T x = var(y t |Y t _ 1 ), T 2 = 
var(y t |X t _!, Y t _i) and T = var(w t |W t _i), where the conditional variance is 
taken to be the variance of the residual about the linear projection which ac- 
counts for the prediction. The measures of linear causality from Y to X, linear 
causality from X to Y, instantaneous linear causality and linear dependence 
were respectively defined to be (Geweke, 1982): 



w t = A(L)e t , 




B(L)wt = £ t , 



(2) 



i 71 Y^x = ln(|S 1 |/|S 2 |), 
i 71 x^Y = ln(|T 1 |/|T 2 |), 
F x . Y = ln(|S 2 |-|T 2 |/|T|), 

F x ,y = IndSxl • l^l/l Y|) = F Y ^ X + + ^x-y- 



(3) 
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The measures of directional linear causality may be decomposed by frequency. 
Let the autoregressive representation for X and Y be: 



B n (L) B 12 (L) 



\B 21 (L) B 22 (L) ] 




(4) 



with B n (0) = I fc , B 22 (0) = l h B 12 (0) = 0, B 21 (0) = 0, var(e lt ) = E 2 , var(e 2t ) = 
T 2 . Eq.(4) is actually a partition form of Eq.(2). Let C = cov(ei t , e 2 t)- Then 
pre-multip lying a transformation matrix 



(5) 




to both sides of Eq.(4), we have the following normalized form: 



/ B 11 (L) B 12 (L) 



\B 2 i(L) B 22 (L) J \y t 



I 




- 


"1 











(6) 



where £i t and £2t are uncorrelated, var(£ 2 t) = T 3 = T 2 — C'S 2 1 C, and 
Bi 2 (0) = but B 2 i(0) 7^ in general. Then Eq.(6) implies the following 
spectral decomposition of the spectral density of X: 



S X (A) = H 11 (A)S 2 HJ 1 (A) + H 12 (A)T 3 H1 2 (A), 



(7) 



where H(A) is the transfer matrix of the normalized autoregressive expression 
in Eq.(6). It is obvious that the spectral density of X is decomposed into 
an intrinsic part and a causal part, so the measure of linear causality was 
suggested as (Geweke, 1982): 



/ Y ^x(A) = ln-~ 



|S X (A)| 



H 11 (A)S 2 HJ 1 (A)| 



(8) 



There is also a convergence relation between the measures in the time and 
frequency domains: 



1 n 



(9) 



In the above, a Granger causality measure between two (or two blocks of) 
time series was given. Before Geweke presented his logarithm version, Pierce 
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(Pierce, 1979) had introdued a R 2 measure which simply takes the ratio of the 
variances of two prediction errors. The value of Pierce's R 2 measure is within 
[0, 1] which is more convenient for comparison with correlation coefficients. 
However, Geweke's logarithm version has better statistical properties. There 
has also been a measure based on autoregressive moving average (ARM A) mod- 
els (Boudjellaba et al, 1992). 

Granger causality analysis has been employed in a number of studies of neural 
data (Bernasconi and Konig, 1999; Bernasconi et al., 2000; Kaminski et al., 
2001; Hesse et al., 2003; Harrison et al, 2003; Brovelli et al., 2004; Roebroeck 
et al., 2005). The procedure described above has potential applications in 
these types of study. For those cases where more than two scalar/block time 
series recordings are available, the procedure may be performed to identify 
further patterns of neural interaction after a more traditional pairwise analy- 
sis. We now consider two simple simulations to illustrate situations in which 
conditional Granger causality analysis is important. 

2.1 Example 1: The case of differentially delayed driving 

We first consider a simple system consisting of three variables, each represent- 
ing an AR process: 



z(t) = /iz(t - 1) + x(t - 2) + e(t). 

where < 1 is a parameter, and £(t), i]{t), e(t) are independent white noise 
processes with zero mean and variances o\ 1 <j\ 1 o~\ 1 respectively. The system 
configuration is illustrated in Fig. 1(a) where x drives y after the delay of one 
time unit and x drives z after the delay of two time units. We note that the 
time unit here is arbitrary and has no physical meaning. To be consistent with 
the data presented later we assume that the sample rate is 200 Hz. In other 
words each time unit is 5 millisecond. 

We performed a simulation of this system, with jj, = 0.5, <7i = 1,<7 2 = 0.2, 
and (T3 = 0.3, to generate a data set of 500 realizations, each 100 points long. 
Then, assuming no knowledge of Eq.(10), we fit multivariate autoregressive 
(MVAR) models (Ding et al., 2000) to the generated data set for each pairwise 
combination of variables x, y, and z, and calculated the frequency-domain 
Granger causality for each pair in each direction, as shown in Fig. 1(b). In the 
top two rows of this figure, we see non-zero Granger causality values across the 
spectra of x — > y and x — > z and zero values across the spectra of y — > x and 



x(t) 

y(t) 



x(t 



1)+ 77(f) 



(10) 
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z — > x. These results are indicative of the true unidirectional causal driving 
of y and z by However, we also see results in the third row of Fig. 1(b) 
which appear to indicate unidirectional causal driving of z by y. In fact, we 
know from the system configuration that this apparent driving is due to the 
common influence of x on both y and z but with different time delays. This 
mistaken identification of an indirect influence as being a direct one suggests 
the need for the conditional Granger causality measure. 

2.2 Example 2: The case of sequential driving 

Next we consider another simple system, again consisting of three AR pro- 
cesses: 



z(t) =fiz(t-l)+y(t-i) + e(t). 

This system configuration consists of sequential driving from x to y, then from 
y to z as illustrated in Fig. 2(a). The same numbers of realizations and data 
points were generated for the same parameter values, and MVAR models were 
again fit to the data. The results of Granger causality analysis in Fig. 2(b) 
show an apparent unidirectional causal driving of z by x that is in fact due 
to the indirect influence through y. Again, the mistaken identification of an 
indirect influence as being direct suggests the need for the conditional Granger 
causality measure. 

Note that, although the systems in the above two examples are very differ- 
ent, the results of pairwise Granger causality analysis seen in Figs. 1(b) and 
2(b) are essentially the same, indicating that the analysis could not distin- 
guish between the two systems. These two examples, although simple, thus 
plainly demonstrate that the pairwise measure of Granger causality by itself 
may be insufficient to reveal true system relations. We now describe the con- 
ditional Granger causality as a potentially useful tool for disambiguating such 
situations. 



3 Geweke's Measure of Conditional Feedback Causality 

Now suppose that w t has been decomposed into three vectors x t , y t and z t 
with k, I and m dimensions, respectively: w t = (xj, y' t , z' t )'. The measure given 



x(t) 

y(t) 



x(t 



1)+ 77(f) 
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by Geweke for the linear dependence of X on Y, conditional on Z, in the time 
domain (geweke, 1984) is: 



Y->X|Z 



In 



var(x f |Xt_i, Z f _i) 
var(x t |X t _i, Y t _i, Z t _i) ' 



(12) 



which is consistent with Granger's definition of a prima facie cause (Granger, 
1980). 

Time series prediction is achieved by the fitting of MVAR models. In order 
to implement Eq.(12), two MVAR models are involved. One is the following 
two-variable MVAR model: 



D n (L) D 12 (L) 



\^D 21 (L) D 22 (L) J 









(**) 











(13) 



with the normalization D n (0) = I, D 22 (0) = I, D 12 (0) = 0, and cov(© t , * t ) = 
imposed in order to yield the frequency decomposition of the conditional 
dependence. The normalization can be achieved by using a transformation 
matrix like Eq.(5). 

The other MVAR model used for deriving the frequency decomposition of the 
conditional dependence is the following three-variable MVAR model: 



Bn(L) 


Bi 2 (L) 


B 13 (^ 




Is \ 


B 2 i(L) 


B 22 (L) 


B 23 (L) 


yt = 


~ £yt 


B 3 i(L) 


B 32 (L) 


B33(L) 


\z t ) 


\e zt ) 



(14) 



with normalization imposed too. The explicit formula of the transformation 
matrix to normalize the MVAR model of three time series is given in the 
Appendix. 

Based on the relations of different variances, Geweke derived the following 
important relation of the conditional causality in the time domain (Geweke, 
1984): 



'y-ociz 



Y*^0- 



(15) 



The same relationship is satisfied in the frequency domain (Geweke, 1984): 
f Y ^x|z(A) = fYv^e(A). (16) 
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In order to get fv*- > e(A), we need to decompose the variance of into the 
frequency domain. To do so, we write Eq.(13) and Eq.(14) in the frequency 
domain: 



^X(A) \ _ I G XX (X) G XZ (X)\ 1 0(A) 
Z(A) ) \G ZX (X) G Z2 (A) j \*{\) 



(17) 



^X(A)^ 



Y(A) 
Z(A) 



H xa; (A) H^A) ~H XZ (\) 
H ya; (A) H W (A) Hyz(X) 





'e x (a) n 




Ey(A) 


J 


v E,(A) y 



If the spectra of X(A) and Z(A) from Eq.(17) remain identical to the spectra 
from Eq.(18), then we can substitute Eq.(17) into Eq.(18) to get the following 
equations: 



/ e(A) N 



Y(A) 
*(A) 



^ G XX (X) G^A)^ ^ 



1/ 

G za; (A) G ZZ (X) 

Qa;a;(A) Qxy(X) Qxz(X) 

Qyx(X) Qyy(X) Q yz (X) 

Q**(A) Q zy (X) Q ZZ (X) 



H^A) H^A) H ;CZ (A) 

Hj, x (A) H W (A) Hy Z (X) 

^~H_ ZX (X) Hzy(X) H 2Z (A) 

' / E X (A) N 



\ 


'v x (x)^ 




E»(A) 


/ 


v E,(A) y 



Ey(A) 
E 2 (A) 



where Q(A) = G 1 (A)H(A). From the first equation of Eqs.(19), the spectrum 
of is decomposed into the following three obvious parts: 

S e (A) = Q XX (X)X XX Q* XX (X) + Q xy (X)X yy Q*(X) + Q^(A)S^Q; z (A).(20) 



(19) 



Therefore the measure of causality from to may be described as: 

|Se(A)| 



fv*^e(A) = In 



Qa;a;(A)S :ca ;Q* ;E (A)| 



(21) 



where Se(A) is actually the variance of ®t, namely Se, since @t is white noise 
in Eq.(13). Considering the relation of Eq.(16), we could get the conditional 
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causality as [11]: 

^ ^ \QxxW^xxQ xx {^)\ ' 

In the above derivations, the assumption that the spectra of X(A) and Z(A) 
coming from Eq.(17) and from Eq.(18) are identical is actually very hard to 
satisfy numerically due to practical estimation errors. As an example of this 
problem, consider Fig. 6, where the dashed curves result from performing 
Geweke's conditional casuality procedure. Note that the negative values seen 
here have no interpretation in terms of causality. (A detailed description of 
Fig. 6 is given in a later section.) In the following section, we introduce the 
partition matrix technique to overcome this problem. 



4 Partition Matrix Improvement 



For three blocks of time series x t , y t , z t , we can fit a three-variable MVAR 
model as in Eq.(14) and we can also derive its frequency domain expression 
as in Eq.(18). From Eq.(18), writing an expression only for X(A) and Z(A) 
(making partitions) we have: 



X(A) 



H ra (A) H xz (\) 



/ E X (A) 



\Z(X)J \U ZX (X) H 22 (A) 



E,(A) 



(23) 



where E X (A) and E 2 (A) have the following moving average expression: 



E*(A) 



X E X (A)^ 



+ 



H xx (\) H xz (\) 



-i 



Hxj/(A) 



^e 2 (a); \v z (\) ) \n zx (\) h 22 (a); \H zy (X) 



Ey(A). (24) 



We realize that Eq.(23) is actually a summation of multiple ARMA processes, 
and that the summation of several ARMA processes is still an ARMA pro- 
cess (Granger and Morris, 1976; Harvey, 1993). However, an unambiguous 
representation of the general multivariate ARMA process for the summation 
is unknown, although a general univariate ARMA model could be obtained 
through a specific procedure (Maravall and Mathis, 1994). Alternatively, we 
adopt the following procedure to evaluate the conditional Granger causality. 

Letting 
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H XJ/ (A) 



H XX (A) H xz (\) 



1 / 

V 



Hxj/(A) 
H 2J/ (A) 



^H 2J/ (A) y yH 2X (A) H 22 (A) 
we get the covariance matrix of the noise terms given in Eq.(24): 



£(A) 



■^-'xx 






^H XJ/ (A) \ 








s 22 y 




v H 2y (A) y 



+ s 



Hxj,(A) 
H zy (A) 



(h;,(a) h^a 

(h:,(A) H^A)). (25) 



This covariance matrix is no longer a real matrix, but it is a Hermite ma- 
trix, i.e. £ X z(A) = S* X (A). Therefore we can use the following transformation 
matrix to normalize the bivariate model of Eq.(23): 



/ 

V 



Ik 

S X2 (A) 







\ 



(26) 



Therefore, in correspondence with the normalized form in Eq.(17), the transfer 
matrix G(A) is now: 



G(A) 



^H 2a ,(A) H 22 (A) 



it o\ 

Sxz(A) 



(27) 



I 



Taking the expansion form of this G(A) matrix to get matrix Q(A) = G _1 (A)H(A), 
and considering S e = E ra , where Y> xx comes from Eq.(25), we can still use 
Eq.(22) to get the conditional causality. 



5 Applications to Simulated and Neural Field Potential Data 



5.1 Application to Simulated Data 



We performed conditional Granger causality analysis on the delay driving 
and sequential driving systems presented above in Section 2. For the delay 
driving case (Section 2.1 and Fig. 1), the Granger causality spectrum from y 
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to z, conditional on x, is presented in Fig. 3. It is obvious from Fig. 3 that the 
conditional Granger causality measure eliminated the indirect causal influence 
of y on z which appeared in Fig. 1(b). For the sequential driving case (Section 
2.2 and Fig. 2), the Granger causality from x to z, conditional on y, is also 
presented in Fig. 3. Clearly, the indirect causal influence from x to z, which was 
indicated in Fig. 2(b), was also eliminated by use of the conditional Granger 
causality. 

In both cases, we have seen that conditional Granger causality analysis elim- 
inated indirect causal influences that inadvertently resulted from application 
of the pairwise Granger causality measure. Knowing the system equations in 
these examples allowed us to verify that the conditional Granger causality 
measure yielded a truer depiction of the system relations. We now consider 
how the conditional Granger causality measure may provide the same benefit 
in the analysis of real neural data. 



5.2 Application to Neural Field Potential Data 

Field potential data were recorded from two macaque monkeys using transcor- 
tical bipolar electrodes at 15 distributed sites in multiple cortical areas of one 
hemisphere (right hemisphere in monkey GE and left hemisphere in monkey 
LU) while the monkeys performed a GO/NO-GO visual pattern discrimination 
task (Bressler et al., 1993). The presence of oscillatory field potential activity 
in the beta (14-30 Hz) frequency range was recently reported in the sensori- 
motor cortex of these monkeys during the prestimulus period (Brovelli et al., 
2004). In that study, Granger causality analysis was performed for all pairwise 
combinations of sensorimotor cortical recording sites. In both monkeys, signif- 
icant Granger causal influences were discovered from primary somatosensory 
cortex to both primary motor cortex and inferior posterior parietal cortex, 
with the latter area also exerting Granger causal influences on primary motor 
cortex. 

In monkey GE, the possibility existed that the causal influence from the pri- 
mary somatosensory (Soma) site to one of the inferior posterior parietal sites 
(in area 7a) was actually mediated by another inferior posterior parietal site 
(in area 7b) (Fig. 4(a)). We therefore used conditional Granger causality anal- 
ysis to test the hypothesis that the Soma — > 7a influence was mediated by the 
7b site. In Fig. 4(b) is presented the pairwise Granger causality spectrum from 
the Soma site to the 7a site (Soma — > 7a, light solid curve), showing signifi- 
cant causal influence in the beta frequency range. Superimposed in Fig. 4(b) is 
the conditional Granger causality spectrum for the same pair, but with the 7b 
site taken into account (Soma — > 7a\7b, dark solid curve). The corresponding 
99% significance thresholds are also presented (light and dark dashed lines 
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which overlap each other). These significance thresholds were determined us- 
ing a permutation procedure that involved creating 500 permutations of the 
field potential data set by random rearrangement of the trial order. Since the 
test was performed separately for each frequency, a correction was necessary 
for the multiple comparisons over the whole range of frequencies. The Bon- 
ferroni correction could not be employed because these multiple comparisons 
were not independent. An alternative strategy was employed following (Blair 
and Karniski, 1993). The Granger causality spectrum was computed for each 
permutation, and then the maximum causality value over the frequency range 
was identified. After 500 permutation steps, a distribution of maximum causal- 
ity values was created. Choosing a p-value at p — 0.01 for this distribution 
gave the thresholds shown in Fig. 4(b), (c) and Fig. 5(b) in dashed lines. 

We see from Fig. 4(b) that the conditional Granger causality is greatly re- 
duced in the beta frequency range and no longer significant, meaning that 
the causal influence from the Soma site to the 7a site is most likely an in- 
direct effect mediated by the 7b site. This conclusion is consistent with the 
known neuroanatomy of the sensorimotor cortex (Felleman and Essen, 1991) 
in which area 7a is connected with area 7b, but not directly with the primary 
somatosensory cortex. 

From Fig. 4(a) we see that the possibility also existed that the causal influence 
from the Soma site to the primary motor (Mot) site in monkey GE was 
mediated by the 7b site. To test this possibility, the Granger causality spectrum 
from Soma to Mot (Soma — > Mot, light solid curve in Fig. 4(c)) was compared 
with the conditional Granger causality spectrum with 7b taken into account 
(Soma — * Mot 1 7b, dark solid curve in Fig. 4(c)). In contrast to Fig. 4(b), we 
see that the beta-frequency conditional Granger causality in Fig. 4(c) is only 
partially reduced, and remains well above the 99% significance level. In Fig. 
5(a), we see that the same possibility existed in monkey LU of the Soma to 
Mot causal influence being mediated by 7b. However, just as in Fig. 4(c), we see 
in Fig. 5(b) that the beta-frequency conditional Granger causality for monkey 
LU is only partially reduced, and remains well above the 99% significance 
level. 

The results from both monkeys thus indicate that the Granger causal influ- 
ence from the primary somatosensory cortex to the primary motor cortex was 
not simply an indirect effect mediated by area 7b. However, we further found 
that area 7b did play a role in mediating the Soma to Mot causal influence 
in both monkeys. This was determined by comparing the means of bootstrap 
resampled distributions of the peak beta Granger causality values from the 
spectra of Soma — > Mot and Soma — > Mot\7b by Student's t-test. The signifi- 
cant reduction of beta-frequency Granger causality when area 7b is taken into 
account (t = 17.2 for GE; t = 18.2 for LU, p «< 0.001 for both), indicates 
that the influence from the primary somatosensory to primary motor area was 
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partially mediated by area 7b. Such an influence is consistent with the known 
neuroanatomy (Felleman and Essen, 1991), which shows direct connections 
between area 7b and both primary motor and primary somatosensory areas. 

As a final demonstration of the value of using the partition matrix method 
outlined in this paper to compute conditional Granger causality, we present 
in Fig. 6 a direct comparison of our improved procedure (solid) with Geweke's 
original procedure (dashed) for the Soma — >■ Mot\7b spectra of monkey GE 
(Fig. 6(a)) and monkey LU (Fig. 6(b)). From much previous experience work- 
ing with this field potential data (Brovelli et al., 2004), we know that spectra 
from these cortical areas typically have a single peak in the beta frequency 
range. Geweke's original method is clearly seen to be deficient in these exam- 
ples not only by the multiple peaks and valleys across the spectra, but also 
by the negative values, which have no physical interpretation. We thus are 
confident that the partition matrix technique is a potentially valuable tool to 
be used in the investigation of conditional Granger causality relations between 
neural signals. 
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Appendix: Transformation matrix to normalize a model of three 
time series 

Since the MVAR model, such as in Eq.(14), is usually not normalized, the noise 
terms could be correlated with each other. Let us assume that the covariance 
matrix is given by: 







Y-'xy 


^12 




s = 


y 

^yx 


y 

^yy 


^yz 






y ^ZX 


^zy 


^zz 


) 



In order to make the first noise term independent, we could use the following 
transform: 
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Ik 








— ^zx^xx 


h 








Pl = 



Then the covariance matrix for the transformed noise terms is: 







'yy 





y y-ly y 

^yx^xx ^xy ^yz 





y y-ly 

^yx^xx ^xz 



\ 



~E zy — Yi zx Ti x lYi xy Yi zz — T, zx T, xx l 'E xz 



Again, to make the second and third noise terms independent, the following 
transformation may be made: 



P 2 = 



Ik 

1,0 

— (^zy ~ ^zx^xx ^xy){^yy ~ ^yx^xx ^xy) Im 



Therefore the whole transformation matrix needed to make all three noise 
terms independent is: 



P = P 2 Pi 
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Figure Caption 



Figure 1: (a) Schematic illustration of a simple delay driving system; (b) 

Granger causality spectra by pairwise analysis. 
Figure 2: (a) Schematic illustration of a simple sequential driving system; 

(b) Granger causality spectra by pairwise analysis. 
Figure 3: Conditional Granger causality spectra from y to z, conditional on 

x, in the delay driving example (Fig. 1), and from x to z, conditional on y, 

in the sequential driving example (Fig. 2), showing that the indirect effects 

indicated in Figs. 1(b) and 2(b) have been eliminated. 
Figure 4: Granger causality analysis of field potential data from monkey GE. 

(a) Schematic illustration of significant Granger causal influences in the beta 
frequency band among sites in sensorimotor cortex, (b) The reduction of the 
conditional Granger causality spectrum (dark solid line) below the signifi- 
cance threshold indicates an indirect influence of the primary somatosensory 
site on the 7a site, (c) The fact that the conditional Granger causality spec- 
trum (dark solid line) does not fall below the significance threshold shows 
that there is a direct influence of the primary somatosensory site on the 
primary motor site. That the conditional spectrum is significantly below 
the pairwise spectrum indicates that there is an additional indirect effect 
mediated by 7b. 

Figure 5: Granger causality analysis of field potential data from monkey 
LU. (a) Schematic illustration of significant Granger causal influences in the 
beta frequency band among sites in sensorimotor cortex, (b) As in Fig. 4(c), 
the conditional Granger causality spectrum (dark solid line) does not fall 
below the significance threshold, indicating that there is a direct influence 
of the primary somatosensory site on the primary motor site. Again, the 
conditional spectrum is significantly below the pairwise spectrum, showing 
that there is an additional indirect effect mediated by 7b. 

Figure 6: Comparison of Geweke's original method for computing condi- 
tional Granger causality spectra (dashed) with our partition matrix proce- 
dure (solid), (a) Soma — > Mot \ 7b in monkey GE, corresponding to Fig. 4(c); 

(b) Soma — > Mot\7b in monkey LU, corresponding to Fig. 5(b). Note that 
for both monkeys, Geweke's original method suffers from having multiple 
peaks and valleys across the spectra (believed to be artifactual), and also 
from having negative values, which have no physical interpretation. 
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